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The  specification  of  the  surface  heat  flux  is  essential 
for  synoptic  and  seasonal  prediction  of  the  upper  ocean 
thermal  structure.  Estimates  of  the  surface  heat  flux  have 
been  prepared  for  the  central  North  Pacific  during  January 
1976  through  April  1979  using  archived  fields  from  the  Fleet 
Numerical  Oceanography  Center  {FNOC)  hemispheric  atmospheric 
prediction  model.  Monthly  accumulations  of  the  surface  heat 
flux  are  compared  with  the  change  in  heat  content  above  230 
m  derived  from  temperature  analyses  of  the  North  Pacific  Ex¬ 
periment  TRANSPAC  ship-of- opportunity  program.  Systematic 
differences  are  found  between  the  accumulated  heat  flux 
fields  and  the  oceanic  heat  content  change.  Some  of  the 
differences  are  due  to  excessively  large  changes  in  ocean 
heat  concents  above  a  fixed  level.  However,  our  earlier 
studies  have  suggested  a  bias  of  excessive  upward  surface 
heat  flux,  especially  along  the  southern  boundary  of  the  do¬ 
main.  Assuming  local  heat  balance  over  a  36-monfh  period,  a 
correction  field  to  the  FNOC  surface  heat  flux  estimates  is 
derived.  Separate  correction  fields  for  the  heating  and 
cooling  seasons  demonstrate  a  seasonal  variation  in  the  ac¬ 
cumulated  heat  flux  versus  heat  concent  change  values. 
Thus,  six  bi-monthly  correction  fields  to  be  added  to  the 
FNOC  heat  fluxes  are  prepared  to  enable  these  heat  fluxes  to 
be  used  for  ocean  prediction. 


1.  IHTB0D0CTI3H 


Prediction  of  the  upper  ocean  the  real  structure  re¬ 
quires  a  specification  of  the  surface  heat  fluxes.  The  time 
scale  on  which  the  sensible  and  latent  fluxes,  plus  the  in¬ 
coming  and  outgoing  radiative  fluxes,  must  be  specified  is 
dependent  on  the  oceanic  phenomena  and  the  type  of  model. 
Heat  flux  values  resolved  on  a  time  scale  of  3  hours  or  less 
must  be  provided  if  the  diurnal  variation  in  the  atmospheric 
forcing  is  an  important  consideration  in  the  evolution  of 
the  upper  ocean  thermal  structure  (Garwood,  1977).  Consid¬ 
erable  evidence  has  been  accumulated  {e.g.,  Elsberry  and 
Camp,  1978;  Camp  and  Elsberry,  197  8;  Elsberry  and  Raney, 
1978;  Elsberry  and  Garwood,  1978;  and  others)  that  the  upper 
ocean  thermal  structure  responds  significantly  to  atmospher¬ 
ic  synoptic-scale  forcing.  On  seasonal  time  scales,  the 
surface  heat  flux  is  accumulated  in  the  upper  ocean  layers 
during  spring  and  summer,  and  subsequently  removed  during 
autumn  and  winter. 

Ocean  mixed  layer  models  provide  a  means  of  demonstrat¬ 
ing  the  two  ways  in  which  the  surface  heat  fluxes  affect  the 
prediction  of  ocean  thermal  structure  profiles  (Niiler  and 
Kraus,  1977).  Consider  a  well-mixed  layer  of  variable 
depth,  h,  with  temperature,  T.  The  heat  content  (H)  per 
unit  area  in  a  flayer  of  depth  (h)  is 

H  *  J  yCp  T  (Z)  dz  (1) 

where  is  the  density,  Cp  is  the  specific  heat,  and  z  is 
the  depth  below  the  surface.  The  change  in  heat  content  of 
the  mixed  layer  is 

<3  (pCpTh)  =  pCfhdT  +  pCfTdh 
dt  d  t  dt  ' 

because  the  changes  in  p  a  nd  Cp  are  relatively  small.  Solv¬ 
ing  for  the  temperature  change  from  the  first  term  on  the 
right  side  of  (2)  gives 

dT  _  1  a  (Th)  T  dh  1  dH 

dt  h  a*  h  dt  '  ?cph  it  *  (3) 

It.  mixed  layer  models  (e.g.,  Garwood,  1977),  the  las t  term 
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The  vertical  current  speed,  W  (-h)  ,  at  the  base  of  the  mixed 
layer  will  be  neglected.  This  relationship  indicates  than 
entrainment  mixing,  and  thus  the  downward  heat  flux  at  base 
of  the  layer,  is  only  associated  with  deepening  layers.  The 
fraction  of  turbulent  kinetic  energy  that  is  available  for 
entrainment  mixing  in  the  Garwood  model  is  dependent  on  both 
the  surface  friction  velocity  and  the  surface  buoyancy  flux, 
which  is  determined  partly  by  the  sue  face  heat  flux.  Conse¬ 
quently,  the  surface  heat  flux  contributes  to  both  terms  in 
(3).  The  effect  in  the  first  term  is  reflected  directly  in 
the  change  ir.  heat  content  due  to  the  surface  flux.  The 
thermal  structure  is  strongly  dependent  on  the  vertical  re¬ 
distribution  of  the  heat  via  entrainment  mixing,  which  is 
partly  due  to  the  upward  surface  heat  flux. 

We  conclude  that  specification  of  the  surface  heat  flux 
is  an  essential  factor  for  prediction  of  anomalous  ocean 
thermal  structure.  In  the  next  section,  we  will  briefly  re¬ 
view  the  methods  available  for  estimating  the  surface  hea* 
flux  over  the  ocean.  In  the  following  section,  we  examine 
the  role  cf  the  surface  heat  flux  in  the  oceanic  heat  budg¬ 
et.  We  then  derive  a  correction  field  to  remove  a  bias  in 
the  surface  heat  flux  that  would  be  detrimental  to  ocean 
thermal  structure  prediction. 


2.  SPECIFICATION  OF  THE  SUB  FACE  HEAT  FLUX 

The  hulk  aerodynamic  method  for  calculating  the  sensi¬ 
ble  and  latent  heat  fluxes  involves  the  surface  wind  field 
and  the  iifferences  in  temperature  and  specific  humidity  be- 
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tween  the  sea  and  the  air.  Calculations  of  tha  incoming  and 
outgoing  radiative  fluxes  require  knowledge  of  the  cloud 
cover  as  well  as  the  solar  altitude.  An  example  of  a  system 
for  estimating  the  surface  heat  flux  from  ship  observations 
is  given  by  Clark,  et  al.  (1  974),  Bunker  (1976)  and  Clark 
(1981).  A  common  approach  is  to  calculate  the  surface  flux 
from  each  ship  report  and  then  average  over  some  space  and 
time  interval.  For  example,  all  ship  observations  taken 
during  a  month  within  5°  latitude  and  longitude  may  be  used 
to  represent  the  average  heat  flux  in  that  domain.  Husby 
(1980)  has  compared  the  anomalous  heat  flux  over  six-month 
periods  at  Ocean  Weather  Station  V  (34*N,  164°E)  with  the 
estimate  based  on  merchant  ship  observations  within  a  4 
quadrangle.  He  concludes  that  the  merchant  ship  reports 
must  be  carefully  screened  prior  to  the  heat  flux  computa¬ 
tions. 

The  accuracy  of  this  method  is  clearly  dependent  on  the 
number  of  observations.  There  is  no  assurance  that  the  re¬ 
ports  will  be  randomly  distributed  in  space  or  time.  One 
problem  appears  to  be  a  "fair  weather  bias"  because  the 
ships  tend  to  avoid  bad  weather.  3arnett  (1981)  indicates 
that  a  set  of  heat  fluxes  in  the  central  North  Pacific  Ocean 
estimated  by  this  method  had  a  bias  of  30-45  W  m"'1  ,  which 
would  result  in  an  excessive  estimation  of  the  heat  flux 
into  the  ocean. 

An  estimate  of  the  monthly  heat  flux  is  not  sufficient 
for  short-term  ocean  prediction,  which  requires  information 
on  diurnal,  or  at  least  synoptic,  time  scales.  Synoptic 
maps  of  heat  flux  are  difficult  to  analyze  because  of  the 
sparsity  of  ship  observations.  An  indirect  method  is  pro¬ 
posed  here.  The  atmospheric  prediction  models  also  require 
a  calculation  of  the  heat  flux  at  the  ocean  surface.  The 
heating  package  for  the  Fleet  Numerical  Oceanography  Center 
(FNOC)  model  involves  bulk  aerodynamic  calculations  of  the 
latent  and  sensible  heat  fluxes,  plus  radiative  fluxes  that 
are  a  function  of  the  model-estimated  cloudiness  (Kesel  and 
Winnirghcff ,  1972)  . 


Our  basic  hypothesis  is  that  the  FNOC  atmospheric  pre¬ 
diction  model  heating  package  can  provide  the  heat  fluxes 
necessary  for  ocean  prediction  experiments.  Gallacher 
(1979)  has  described  the  method  used  for  extracting  hourly 
heat  flux  estimates  from  the  FNOC  archives.  Elsberry,  Gal- 
lacher  and  Garwood  (1979)  have  used  these  heat  fluxes  to 
predict  the  ocean  thermal  structure  changes  during  the  au¬ 
tumn  of  1976.  Budd  (1980)  also  used  these  fields  in  an  at¬ 
tempt  to  predict  the  spring  transition  from  the  winter  to 
the  summer  regimes  in  the  central  North  Pacific  Ocean.  Budd 
found  a  systematic  bias  with  a  too  large  upward  heat  flux 
near  30#N.  Recently,  Steiner  (1581)  also  found  a  systematic 
bias  in  these  heat  flux  estimates  in  the  region  between  Ha¬ 
waii  and  San  Francisco. 

The  purpose  of  this  paper  is  to  describe  a  correction 
field  to  the  FNOC  surface  heat  flux  fields  to  permit  their 
use  in  ocean  prediction  experiments.  As  in  the  case  of  Budd 
(1980)  and  Steiner  (1981),  the  approach  is  to  derive  a  cor¬ 
rection  field  that  assures  a  long  term  (annual  or  longer) 
heat  balance  in  the  upper  ccean.  Analyses  by  White  and  Bern¬ 
stein  (1979)  of  the  ocean  thermal  structure  observations  in 
the  TRANSPAC  ship-of-opportunity  program  are  used  to  calcu¬ 
late  the  time  changes  in  oceanic  heat  content.  These  are 
compared  with  the  accumulated  surface  heat  fluxes  derived 
from  the  FNOC  files.  We  then  derive  the  correction  to  the 
monthly  surface  heat  flux  that  is  necessary  to  assure  local 
heat  balance  over  the  1976-1978  period.  In  the  next  section 
we  review  briefly  some  tests  of  the  local  heat  balance  as¬ 
sumption.  It  is  only  on  the  very  long  time  interval  that  we 
assume  local  heat  balance.  We  do  not  require  or  insure  lo¬ 
cal  heat  balance  cn  monthly  or  shorter  time  intervals  over 
which  we  are  doing  prediction  experiments. 

3.  OCEAN  HEAT  BUDGET  STUDIES 

The  general  purpose  in  This  section  is  no  determine  one 
relative  impornar.ee  of  nhe  surface  heat  flux  in  the  upper 
ocean  he  an  budgen.  Our  primary  inneresn  is  in  coen-ccear. 
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regimes  that  would  be  similar  to  the  conditions  in  the  North 
Pacific  Experiment  (NORPAI)  Anomaly  Dynamics  Study  (ADS)  re¬ 
gion.  We  do  not  consider  regions  of  boundary  currents  or 
near-equatorial  areas  in  which  horizontal  and  vertical  ad- 
vecticn  are  likely  to  be  significant.  The  space  and  time 
scales  of  interest  are  greater  than  1  000  km  and  one  month. 

Many  comparisons  have  been  made  between  the  seasonal 
changes  of  heat  content  of  the  upper  ocean  and  the  accumu¬ 
lated  heat  flux  at  the  surface.  Bryan  and  Schroeder  (1960) 
compared  the  heat  content  calculated  from  North  Atlantic  BT 
data  with  the  surface  heating  estimated  by  Budyko  (1955). 
They  found  that  the  surface  heating  on  a  seasonal  basis  was 
about  20$  less  than  the  change  in  heat  content  in  the  region 
between  20®N  and  about  50*N.  By  contrast,  Bathen  1971) 
found  that  the  surface  heat  exchange  estimates  of  rtki 
(1966)  could  account  for  only  29$  of  the  iocal  a  hly 
change  in  heat  content  in  the  North  Pacific  Ocean.  Gi  and 
Niiler  (  1973)  propose  that  the  heat  input  averages  er 

large  areas  and  times  is  mainly  stored  locally,  and  horizon¬ 
tal  advection  by  the  mean  flow  is  not  particularly  impor¬ 
tant.  They  also  cite  comparisons  by  Tabata  (1955)  and  Rob¬ 
inson  (1966)  using  Ocean  Weather  Ship  data  which  suggest 
that  most  of  the  heat  input  changes  are  stored  locally. 
Gill  and  Niiler  further  suggest  that  the  inaccuracy  of  the 
heat  flux  estimates  over  large  areas  away  from  Ocean  weather 
Ships  is  the  likely  cause  cf  some  of  the  departure  from  lo¬ 
cal  heat  balance. 

A  careful  study  of  the  upper  ocean  heat  budget  near  OWS 
?  based  on  two  weeks  of  high-quality  observations  during  the 
Mixed  Layer  Experiment  (MILS)  has  been  reported  by  Davis,  ec 
ai.  (198  1).  They  found  that  a  one- dimensional  upper  layer 
heat  budqet  may  be  closed  acceptably  if  the  temperature  is 
well-sampled,  and  if  the  vertical  velocity  in  the  seasonal 
thermocline  is  also  taken  into  account. 

A  recent  study  by  Barnett  ( 1 9 3 1  •  used  AXBT's  along  153 

and  1 7  0  *  W  between  30*  and  50*  N  at  approximately  monthly  in¬ 


tervals  tc  estimate  the  heat  budqet.  Barnett  concludes  that 


approximately  90-95%  of  the  variance  in.  the  seasonal  change 
of  heat  storage  can  be  accounted  for  by  the  air/sea  heat  ex¬ 
change.  Horizontal  and  vertical  advectior.  were  of  limited 
and  no  use,  respectively,  in  reproducing  the  seasonal  cycle 
variance.  Barnett  also  concludes  that  prior  heat-budget 
studies  that  used  an  arbitrary  lower  depth  in  the  estimation 
of  heat  content  may  have  inadvertently  included  variance  on 
different  space  and  time  scales.  Finally,  Barnett  suggests 
that  the  inability  to  explain  adequately  the  anomalous  heat 
content  changes  in  the  central  Pacific  is  likely  due  to  im¬ 
precise  knowledge  cf  the  source  terns  in  the  heat  budget. 

Our  earlier  studies  (Elsberry  efc  al.  ,  1979;  Budd,  1989; 
Steiner,  1981)  have  used  an  oceanic  heat  budget  to  determine 
the  feasibility  of  using  the  FNOC  heat  flux  estimates  for 
ocean  prediction.  In  each  case,  the  imbalances  appear  to  be 
systematic  in  space  and  time.  Large,  short-term  imbalances 
are  likely  due  to  ocean  observational  sampling  errors  and  to 
using  an  arbitrary  lower  depth  in  the  estimate  of  heat  con¬ 
tent  (Barnett,  1981).  However,  the  long-term,  systematic 
error  is  likely  due  to  a  bias  in  the  FNOC  heat  flux  esti¬ 
mates.  In  the  following  sections,  we  derive  an  appropriate 
correction  field  to  be  added  to  the  P  NOC  heat  flux  estimates 
to  remove  the  long  term  bias. 

4.  CALCULATION  CF  THE  OCEAN  HEAT  BUDGET 

Monthly  mean  temperatures  during  1976-1979  have  been 
objectively  analyzed  by  White  and  Bernstein  (  1979)  cased  on 
TRAHSPAC  ship-of-opportur.ity  XBT's  only.  Analyzed  tempera¬ 
tures  are  available  on  a  2°  latitude  by  5 *  longitude  grid  at 
0,  20,  40,  60,  90,  120,  1  5  0,  200,  259  ,  300  and  4C0  a.  The 
general  domain  cf  the  heat  budget  calculations  is  from  30" 
tc  50* N  and  from  170*E  to  130*W.  However,  taa  number  of  ob¬ 
servations  in  the  northwest  and  southwest  corners  of  the 
grid  are  insufficient  to  provide  reliable  estimates  -during 
seme  months.  Consequently,  these  areas  are  eliminated  from 
the  following  analysis,  and  will  appear  as  cross-hat ched 
areas  on  all  figures. 
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To  calculate  the  heat  content  as  in  (1),  the  objective¬ 
ly  analyzed  temperatures  are  first  interpolated  to  5  a 
depths.  The  heat  content  is  actually  calculated  relative  to 
the  200  m  temperature,  with  the  intent  of  removing  some  of 
the  change  ia  heat  content  due  to  vertical  displacements  of 
■the  thermccline.  Barnett  (1981)  has  suggested  that  selec¬ 
tion  of  an  arbitrary  depth  such  as  20  0  m  may  introduce  addi¬ 
tional  variaace  into  the  heat  content  calculation  due  to  in¬ 
ternal  waves,  tides  and  so  forth.  Inspection  of 

month-to- month  heat  content  values  does  reveal  time  changes 
which  are  much  larger  than  can  be  reasonably  accounted  for 
by  surface  heat  fluxes.  Large  positive  time  changes  tend  to 
be  followed  by  large  negative  changes.  Thus,  time  intervals 
of  either  two  or  three  months  have  been  tested  in  the  esti¬ 
mate  of  the  oceanic  heat  content  change, 

=  H(t  ♦  A*)  -  H  (t)  .  (5) 

He  expect  AH  to  be  positive  during  the  period  of  net  down¬ 
ward  surface  heat  flux  (roughly  between  April  and  September) 
and  negative  during  the  period  of  net  upward  flux.  The 
monthly  TEANSPAC  analyses  are  assumed  to  apply  on  the  15th 
of  each  month.  If  At  is  two  months  in  (5)  ,  then  /\H  will 

also  apply  at  the  15th  of  the  month.  For  example,  the  dif¬ 
ference  between  the  March  1976  and  the  January  1976  heat 
content  fields  would  be  applied  at  the  middle  of  February. 
Notice  that  this  is  the  first  possible  value  that  can  be 
calculated  because  the  analyses  begin  in  January  1976.  How¬ 
ever,  if  At  three  months  in  (5),  then  Ad  will  apply  at 
the  mid-pcint  of  this  interval.  Thus,  the  first  possible 
difference  would  be  between  April  1976  and  January  1976,  and 
i t  would  apply  on  1  March  1976.  The  advantage  of  the  longer 
rime  interval  is  that  shorter  period  fluctuations  in  heat 
content  tend  to  be  averaged  out  and  one  obtains  a  mere  con¬ 
sistent  measure  of  the  seasonal  variation.  The  disadvantage 
of  the  longer  time  interval  is  the  greater  inaccuracy  in  the 
finite  difference  approximation  to  a  derivative.  Finally, 
it  is  necessary  tc  extend  the  AH  calculation  into  197?  to 
obtain  three  complete  annual  cycles.  In  the  case  of  A*  -  - 


(3)  months,  the  analyses  through  February  (March)  1979  are 


used. 


Following  Emery  (1976), 
may  be  written  as 

—  -  -  ?cpot  > 

dt  ~  D 


the  basic  heat  balance  equation 

-  V  •  ♦  Q  ,  (6) 


where  H  is  defined  as  in  (1)  with  -h ,  the  mixed  layer 
depth,  replaced  by  D  which  is  a  (constant)  depth  below  the 
main  theraocline.  In  (6)  ,  the  first  term  on  the  right  rep¬ 
resents  the  divergence  of  heat  due  to  horizontal  divergence, 

Wp  H/D,  and  vertical  advection  WpfCpT  .  The  second  term  is 
the  horizontal  advection  of  heat  and  the  last  term  is  the 
surface  heat  flux.  The  suggestion  by  Barnett  (1981)  and 
Davis  et .  al.  (1981)  that  D  be  the  depth  of  an  isotherm  be¬ 
low  the  main  theraocline  rather  than  a  constant  depth  re¬ 
sults  in  the  divergence  term  in  (6)  being  zero,  and  thus 
eliminates  a  source  of  uncertainty  from  the  heat  budget. 
Due  to  the  large  vertical  sampling  interval  in  the  TRANSPAC 
analyses,  we  will  use  D  =  200m  throughout  this  study.  Thus, 
we  assume  the  following  budget  equation 

=  /  QSt  ♦  Residual  *  QT  ♦  Residual  .  (7) 

The  integral  sign  indicates  that  the  air-sea  flux  is  summed 
over  the  same  time  interval  that  is  evaluated.  QT  is 

defined  as  the  surface  heating.  In  the  =  2  month  exam¬ 
ple  above,  the  sum  would  be  from  15  January  through  15 
March.  The  residual  term  in  (7)  includes  the  remaining 
terms  (physical  effects)  in  (6)  that  can  not  be  evaluated, 
plus  the  errors  in  estimating  the  heat  content  changes  and 
the  surface  fluxes.  If  these  observational  errors  are  ran¬ 
dom  over  a  sufficiently  long  interval,  their  contribution 
should  be  averaged  cut.  In  particular,  we  sum  ever  36 
months  (or  a  small  subset  if  a  seasonally  varying  correction 
factor  is  derived)  and  derive  an  averaged  correction  field 
to  be  added  to  the  surface  heat  flux.  However,  t he  physical 
terms  may  not  sum  to  zero.  Persistent  vertical  ana  horizon¬ 
tal  advection  or  diffusive  effects  that  accumulate  during 
the  period  are  thus  also  incorporated  into  ‘•he  correction 


field.  The  inpact  of  assuning  a  local  (one-dimensional) 
heat  balance  via  this  long-tern  correction  field  will  be  as¬ 
sessed  in  the  prediction  experiments. 

The  total  heat  flax  Q  in  (7)  may  be  expressed  as 

C  =  Q*  -  <Qb  ♦  ♦  Q1)  r  ( 
where  the  subscripts  s,  b,  1  and  h  refer  to  solar,  back,  la¬ 
tent  and  sensible  heat  fluxes  through  the  sea  surface.  A 
description  of  the  calculation  of  each  of  these  heat  fluxes 
in  the  FNOC  atmospheric  prediction  model  is  given  in  Kesel 
and  Winninghoff  (1972).  The  values  cf  these  terms  after  a 
one-hour  integration  of  the  model  are  taken  to  apply  at  the 
synoptic  map  time  (00  or  12  GMT).  After  spatial  interpola¬ 
tion  to  the  grid  points  of  the  ocean  analysis,  these  instan¬ 
taneous  values  are  interpolated  to  hourly  values  using  tech¬ 
niques  described  by  Gallacher  (1979).  The  representation  of 
the  solar  flux  is  particularly  involved.  It  is  this  time 
series  of  hourly  values  of  Q4  and  (=  *•  3K  ♦  Qx)  that 
are  accumulated  in  (7).  Simply  summing  the  instantaneous 
12-h  FNOC  values  would  not  produce  the  same  sum  because  we 
have  reconstructed  the  hourly  variation  in  Q. 

An  example  of  a  field  of  near  the  beginning  of  the 
ocean  heating  season  is  given  in  Fig.  1  .  Gver  most  (all) 
of  the  domain,  the  heat  content  change  is  positive,  as  ex- 

(f  -2 

pected.  A  value  of  2  x  10  cal  cm  over  this  two-month 
period  corresponds  to  an  average  increase  in  temperature  of 
2°  C  over  a  depth  of  100  m,  or  4°C  over  a  depth  of  50  m. 
Less  confidence  should  be  given  to  the  values  near  the 
northern  and  southern  boundaries  of  the  domain.  The  primary 
ship-of- opportunity  tracks  lie  on  the  great,  circle  route  be¬ 
tween  Japan  and  the  USA  west  coast  ports  (White  and  Bern¬ 
stein,  1979).  During  the  winter,  the  tracks  tend  to  be  dis¬ 
placed  southward,  which  makes  the  northern  region  sampling 
rather  poor.  Somewhat  the  opposite  effect  occurs  during  the 
summer  as  the  ship  tracks  are  extended  northward.  The  ob¬ 
jective  analysis  procedure  (White  and  Bernstein,  1979)  used 
through  1978  would  generate  values  at  all  gridpoints  regard¬ 
less  cf  whether  any  observations  fell  in  the  vicinitv  of  the 


point.  Because  the  deficient  points  tend  to  lie  toward  the 
north  and  the  south,  the  fictitious  values  tend  to  resemble 
northward  or  southward  extrapolations  from  the  central  re¬ 
gion. 

The  corresponding  integrated  surface  heat  flux  for  this 
period  derived  from  the  FNOC  calculations  is  shown  in  Fig. 
1c.  In  contrast  to  the  fields,  no  extrapolation  of  the 
fields  toward  the  domain  boundaries  is  involved  in  the  QT 
field.  In  the  band  from  40-50°  N ,  the  surface  exchange  is 
approximately  equal  to  the  observed  ocean  heat  content 
change,  which  suggests  an  approximate  local  heat  balance. 
However,  the  net  FNOC  heat  flux  continues  to  be  upward  be¬ 
tween  30-35°  N  during  this  period.  One  can  not  tell  from 
this  diagram  whether  the  daily  values  of  Q3  in  (8)  are  too 
small,  or  if  the  upward  surface  heat  flux  (Qa)  is  too  large. 
The  fact  that  the  maximum  values  are  found  near  the  longi¬ 
tude  for  which  the  boundary  of  the  FNOC  grid  is  tangent  to 
the  eguator  (therefore  closest  to  the  study  region)  suggests 
that  the  problem  may  lie  in  the  boundary  conditions  that  are 
applied. 

The  difference  between  and  the  cumulative  heat  flux 
(Qt)  is  shown  in  Fig.  1  .  A  licjht  filter,  with  weights 
2-4-2  in  latitude  and  1-4-1  in  longitude,  has  been  passed 
over  this  field  to  reduce  small  scale  noise.  Although  there 
are  considerable  areas  with  the  expected  zero  values,  there 
are  other  araas  with  positive  and  negative  values.  In  par¬ 
ticular,  the  residual  term  in  (8)  is  very  large  along  the 
southern  boundary  because  of  the  lack  of  a  net  downward  sur¬ 
face  heat  flux  (Fig.  1  )  . 

The  complete  set  of  monthly  ^3,  2T  and  values 

is  contained  in  the  Appendix.  Since  our  desire  is  to  obtain 
a  correction  field,  if  is  not  appropriate  to  discuss  here 
each  of  these  sets.  One  general  feature  is  that  the  QT 
fields  have  less  east-west  variation  and  more  north-south 
variation  than  do  the  fields.  There  are  several  possi¬ 
ble  causes  for  this  feature:  (a)  the  FNOC  heat  fluxes  may 

net  contain  an  adequate  representation  cf  the  east-west  var- 
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iations  in  the  anomalous  surface  fluxes;  (b)  vertical  adjec¬ 
tive  processes  associated  with  ocean  eddies  are  contained  in 
the  heat  content  values  relative  to  a  fixed  depth;  and  (c) 
errors  due  to  inadequate  observations  or  analysis  of  the 
ocean  heat  content.  The  spring  transition  regime  example  in 
Fig.  1  illustrates  one  of  the  worst  imbalances  in  the  time 
series.  An  excessive  upward  heat  flux  also  occurs  along  the 
southern  boundary  during  autumn.  However,  the  field  is 
negative  (upper  ocean  heat  content  diminishing  in  time)  dur¬ 
ing  this  period.  Consequently,  tha  magnitude  and  sign  of 
the  imbalance  may  have  a  seasonal  dependence  (discussed  fur¬ 
ther  in  the  next  section)  . 

Similar  sets  of  QT  and  ^a-QT  were  also  prepared 
with  A-  =  3  months.  These  fields  (not  shown!  are  not  di¬ 
rectly  comparable  with  the  two-month  set  because  of  the 
half-month  displacement  of  the  central  point.  However,  the 
time  evolution  of  the  ma  j  cr  features  using  ^t  =  3  months 
can  be  easily  associated  with  the  evolution  shown  by  the  At 
=  2  month  series.  Therefore,  only  the  correction  fields 
based  on  the  2-month  differences  will  be  discussed  below. 

5.  INT ERPRETATIOH  OF  THE  HEAT  BUDGET  IMBALANCE 

Perhaps  the  most  important  factor  to  emphasize  regard¬ 
ing  these  calculations  of  the  local  heat  budget  (8)  is  that 
each  of  the  two  terras  is  subject  to  large  errors.  As  noted 
above,  calculating  the  heat  content  changes  relative  to  a 
fixed  depth  aay  produce  oscillations  equivalent  to  a  1 -2*C 
change  averaged  over  the  entire  depth.  A  faulty  XBT  could 
also  be  the  cause  of  a  temperature  bias  with  depth.  How¬ 
ever,  we  would  expect  that  the  objective  analysis  technique 
would  tend  to  eliminate  such  a  bias  if  other  correct  pro¬ 
files  are  in  that  region.  Based  on  our  earlier  studies 
(Elsberry,  et  al. ,  1979  ;  Budd,  1983;  Steiner,  1981),  we  ex¬ 
pect  that  the  monthly  surface  heating  dees  not  have  a  suffi¬ 
ciently  large  seasonal  amplitude.  There  is  also  a  persis¬ 
tent  bias  toward  excessive  heat  loss  to  the  atmosphere  along 
the  southern  boundary  of  the  domain. 


Because  the  heat  budget  disbalances  may  arise  from  ei¬ 
ther  term.  Table  1  was  prepared  to  indicate  the  various  pos¬ 
sibilities  that  may  arise.  It  is  useful  to  separate  the  AH 
into  separate  periods  when  the  ocean  heat  content  is  in¬ 
creasing  cr  decreasing.  Apart  from  small  calculation  er¬ 
rors,  one  would  hope  that  the  corresponding  2T  would  be  po¬ 
sitive  {surface  heating)  or  negative,  respectively,  so  that 
Cases  A-B  or  E-F  would  apply.  One  expects  approximately  an 
egual  distribution  between  A  and  3  or  between  E  and  F  if  the 

errors  are  random.  The  percentage  of  gridpoints  with  month- 
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ly  Ah“Qt  differences  exceeding  0.5  x  10  cal  cm  are  nor 
evenly  distributed  between  cases  A  and  B  or  between  E  and  F 
when  the  uncorrected  QT  is  used  in  the  differences.  There 
is  clearly  a  bias  toward  cases  A  and  E,  which  could  be  at¬ 
tributed  to  excessive  upward  heat  flux  during  both  the  ocean 
warming  and  cooling  periods.  However  an  approximately  equal 
distribution  is  obtained  when  the  corrected  (see  description 
of  six  bi-monthly  correction  fields  in  the  next  section) 
surface  heat  flux  is  used  in  the  differences. 

4  -2 

A  value  of  0.5  x  10  cal  cm  corresponds  to  a  heat  flux 
bias  of  7.0  cal  cm  hr  or  a  temperature  bias  of  .25  C  in  a 
200  m  water  column.  Both  of  these  values  are  within  the  ex¬ 
pected  range  of  instrument  error  for  the  measurements  which 
were  used  in  the  analyses.  Whereas  65&  of  all  the  differ¬ 
ences  exceed  this  criterion  for  the  uncorrected  heat  flux, 
only  42^  exceed  the  criterion  for  the  corrected  heat  flux. 

Cases  C  and  D  in  Table  1  are  labelled  as  drastic  imba¬ 
lances  because  the  A H  an<3.  Qr  are  of  opposite  tendencies.  A 
physical  explanation  for  such  an  event  might  be  the  advec- 
tior.  of  a  warm  (cold)  ocean  eddy  into  the  region  that  nas 
upward  (downward)  surface  heat  flux.  An  example  of  a  compu¬ 
tational  explanation  is  found  along  the  southern  boundary  in 
Fig.  1  where  Q T  <  0  when  A 3  >  3.  The  percentages  in  Table 
1  indicate  -hat  drastic  imbalances  are  relatively  rare  dur¬ 
ing  ocean  cooling  periods  (Case  D  -  1  ^)  when  the  uncorrected 
Qt  is  used  in  the  differences.  However,  this  is  not  the 
case  during  -.he  ocean  warming  periods,  when  a  Large  fraction 


-12- 


of  the  points  have  significant  upward  heat  fluxes  rather 
than  downward.  Since  the  heat  budget  for  the  central  Pacif¬ 
ic  is  in  an  approximate  local  balance  for  seasonal  time 
scales  (Gill  and  Niiler,  1973),  the  percentages  for  Cases  C 
and  D  should  be  less  than  those  for  Cases  A,  B,  E  and  ?. 
This  distribution  is  obtained  when  the  corrected  QT  values 
are  used,  but  not  with  the  uncorrected  2T.  This  is  a  fur¬ 
ther  demonstration  of  the  need  to  adjust  the  FNOC  surface 
heat  fluxes. 

The  alternative  of  a  systematic  bias  due  to  the  AH 
calculation  does  net  appear  likely.  Although  there  are  more 
points  with  AH  positive  (59$)  than  negative  (41$),  this  as¬ 
ymmetry  is  probably  within  the  limits  of  roughly  offsetting 
periods  of  ocean  warming  versus  cooling.  The  imbalances  be¬ 
tween  seasons  in  Table  1  are  not  consistent  with  horizontal 
advection  being  a  primary  cause.  One  would  expect  stronger 
Ekman  advection  effects  during  the  winter  (cooling  season), 
whereas  the  larger  imbalances  are  found  during  tie  summer. 

TABLE  1 
A  H  versus  QT 


mm 

Qt  >  0 

(Downward  surface 
heat  flux) 

Qr  .  <  0  . 

(Upward  surface 
heat  flux) 

„  AH  >  0. 
Ocean  warming 

Case  A  AH-QT  >0 

Inadequate  downward  Q 
Excessive  ocean  warming 
22.4%  (26.8%) 

Case  C  AH-0t  >  0 

Drastic  imbalance 

33.4%  (7.3%) 

Case  B  AH-Qt  <  0 

Excessive  downward  Q 
Inadequate  ocean  warming 
2.9%  (20.0%) 

AkH  <  0 

Ocean  coding 

Case  D  .  AH-Qt  <  0 

Drastic  imbalance 

1.0%  (6.5%) 

Case  E  AH-Qr  >  0 

Excessive  upward  Q 
Inadequate  ocean  coclirc 
3u.6%  (16.7%) 

Case  ?  AH-Qt  <  0 

Inadequate  upward  Q 
excessive  ocean  coolinc 
5.4%  (22.6%) 

1 

% 1 
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Possible  physical  association  for  Dositive  or  negative  heat  budeet 
imbalances  civen  that  the  heat  content  or  integrated  surface  fluxes 
are  oositive  or  neaative.  Percentage  of  the  gridooir.ts  for  the  36 
monthly  maos  that  fall  in  each  cat  ego  rv  are  indicated.  The  oer  cent  aces 
in  parentheses  were  comouted  using  the*  corrected  total  heat' flux. 

Only  points  with  IAH~Qt1  ^  0.5  x'10*  cal  cm“a  for  the  two  month 
difrerence  are  included  m  tae  calculation  of  the  oercert. 
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S.  CALCULATION  OF  COBBBCTIOH  FIELD 


Our  objective  is  to  determine  a  correction  field  that 
may  be  applied  to  the  Q  values  in  (3)  during  the  ocean  pre¬ 
diction  experiments.  Because  the  time  step  in  the  Garwood 
model  is  typically  1  n,  it  is  desirable  that  the  correction 
be  in  units  of  cal  cm”2  hour”1.  This  is  dona  by  converting 
each  of  ^H-QT  fields  to  an  hourly  rate  using  the  proper 
number  of  hours  in  the  time  interval.  The  sign  is  also 
changed  because  the  correction  is  to  be  added  to  the  Q,. 
field. 

It  is  not  within  the  scope  of  this  paper  to  determine 
the  fraction  of  this  correction  which  will  be  applied  to  Qs 
versus  the  ramaining  three  terms  in  (8).  It  might  be  noted 
that  Qs  occurs  only  during  the  daytime  and  this  flux  is  to 
be  distributed  exponentially  with  depth  by  the  prediction 
model.  In  contrast,  the  remaining  three  terms  in  (8)  apply 
only  at  the  surface  and  tend  to  be  positive  (upward) 
throughout  the  day.  Thus,  one  may  expect  considerable  dif¬ 
ferences  in  the  upper  ocean  predictions  as  the  fraction  of 
the  correction  that  is  allcted  to  Qs  is  increased. 

For  simplicity,  one  desires  a  single  correction  field 
as  in  Fig.  2.  The  basic  features  in  this  correction  field 
are  generally  consistent  with  the  pattern  shown  in  Fig.  la. 
In  particular,  the  correction  reduces  the  surface  heat  flux 
(upward  is  defined  to  be  positive)  along  the  southern  bound¬ 
ary.  It  can  be  shown  that  this  correction  is  equivalent  to 
calculating  the  residual  in  (7)  using  the  difference  in  H 
between  Januar y-Februar y  1979  and  Jan uary-February  1976,  and 
the  Qt  for  the  entire  period.  We  may  also  regard  the  single 
correction  field  as  the  adjustment  necessary  during  each 
hourly  time  step  of  an  integration  from  15  January  1976  to 
15  January  1979.  If  the  correction  is  applied  in  this  way, 
we  can  be  assured  of  conservation  of  heat  at  each  gridpoint 
for  the  entire  three  year  period  (assuming  no  vertical  dif¬ 
fusion  and  no  rcund-cff  errors). 

Each  of  the  36  maps  of  ^H-Qr  is  different  from  that 
implied  by  ‘he  single  correction  field.  '«' e  noted  above  that 
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the  seasonal  variation  in  or  Qr  light  be  associated  with 
a  Modulation  in  the  difference  field.  A  separate  correction 
field  is  derived  for  November  through  April  (Pig.  3a)  and 
for  May  through  October  (Fig.  3b).  Each  of  these  correc¬ 
tions  is  summed  over  the  three  years  and  six  months,  so  that 
a  total  of  18  values  is  included  in  a  ach  average.  There  is 
clearly  a  seasonal  variation  in  these  two  correction  fields. 
The  pattern  during  the  heating  season  (Fig.  3b)  is  similar 
to  the  single  correction  field  (Fig.  2),  but  the  values  are 
larqer.  The  pattern  during  the  cooling  season  (Fig.  3a)  is 
quite  different.  Not  only  are  the  values  considerably 
smaller  than  in  Fig.  3b,  there  is  mere  of  a  north-south  ori¬ 
entation  of  the  isclines  in  the  north-central  region. 

Based  on  the  seasonality  in  Fig.  3,  we  also  examined 
further  subdivisions  into  quarterly  or  bi-monthly  correction 
fields.  A  decision  was  then  made  to  adopt  the  six  correc¬ 
tion  fields  in  Fig.  4  a-f.  It  is  felt  that  six  maps  will 
give  a  better  representation  of  the  seasonal  variation.  In¬ 
spection  of  Fig.  4  b-e  reveals  a  basically  east-west  pat¬ 
tern,  whereas  the  remaining  two  maps  exhibit  the  north-south 
orientation.  Even  though  the  seasonal  variation  does  not 
appear  to  be  a  sine  wave,  the  six  maps  appear  to  provide  a 
relatively  smooth  transition  between  the  two  basic  patterns 
in  Fio.  3.  There  are  some  non-sea sonal  features  in  these 
correction  fields.  The  largest  is  found  along  the  western 
side  of  the  domain.  Larger  corrections  are  required  during 
February- Marc h  (Fig.  4a)  and  during  August-Saptember  (Fig. 
4 d) .  The  transition  between  the  correction  field  for  Decem- 
ber-January  (Fig.  4e)  and  for  Febr l  ary- March  is  especially 
noteworthy  along  *he  western  boundary.  The  regions  of  posi¬ 
tive  corrections  (additional  upward  heat  flux  required)  also 
tend  to  be  somewhat  erratic.  However,  these  positive  cor¬ 
rections  are  always  small,  so  they  will  have  little  effect 
on  the  ocean  predictions. 

An  example  of  the  application  of  the  heat  flux  correc¬ 
tion  field  at  a  specific  location  is  given  in  Fig.  5.  The 
seasonal  changes  in  heat  content  calculated  from  the 
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TRANSPAC  analyses  have  relatively  constant  values  during  the 
heating  season.  By  contrast,  the  values  during  the 

cooling  season  tend  to  reach  rather  well-defined  peaks.  The 
surface  heating  agrees  rather  well  with  the  maximum  ^  H  val¬ 
ues  during  each  cooling  season.  However,  the  uncorrected 
heat  flux  during  the  warming  season  is  clearly  deficient, 
especially  during  1978.  Applying  the  heat  flux  corrections 
for  this  location  from  Fig.  4  improves  the  agreement  between 
the  surface  heating  and  the  ^  H.  The  major  feature  to  be 
noticed  is  that  the  (^  H  -  QT)  differences  ace  indeed  sys¬ 
tematic,  so  that  a  single  correction  in  each  two-month  peri¬ 
od  tends  to  improve  all  three  years.  One  finds  periods  in 
which  the  -  QT  remains  large  after  the  correction  has 

been  applied.  This  indicates  that  although  we  have  had  to 
correct  the  surface  heat  flux  we  have  not  forced  the  heat 
budget  to  be  one-dimensional  for  periods  shorter  than  36 
months.  Some  examples  are  the  differences  during  April  - 
June  1976  and  during  July  -  August  1978.  Another  feature  in 
Fig.  5  is  that  the  area  between  the  H  and  the  corrected  qt 
curves  must  sum  to  the  difference  between  the  final  and  ini¬ 
tial  H  values  (February  1979  and  January  1976,  respective¬ 
ly).  This  requirement  is  a  consequence  of  the  local  heat 
budget  assumption  over  the  36-month  period. 

7.  SUMMARY 

We  have  prepared  six  bi-monthly  correction  fields  to  be 
applied  to  the  FNOC  heat  flux  values  to  be  used  for  ocean 
prediction.  The  largest  corrections  are  found  generally  be¬ 
tween  30  and  33  N  during  April  through  November.  The  heat 
flux  bias  is  evidently  not  serious  for  the  atmospheric  pre¬ 
dictions  because  they  are  limited  to  72  h .  However,  such  a 
bias  can  be  disastrous  <3udd,  19901  for  ocean  prediction 
over  monthly  time  periods.  The  correction  fields  are  aver¬ 
aged  over  three  years  (19  76-  1  979)  .  Closure  of  the  local 
heat  tudaet  over  the  entire  period  is  insured  by  the  use  of 
this  correction  field.  However,  this  is  not  true  over  ary 
shorter  time  intervals.  It  particular,  fluctuations  ir.  heat 
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content  m  the  monthly  TRANSPAC  analyses  are  not  accounted 
for  by  the  long-term  correction  field. 

The  success  of  these  correction  fields  can  only  be 
judged  by  thair  application  in  the  ocean  prediction  experi¬ 
ments.  Inclusion  of  the  corrections  should  result  in  im¬ 
proved  predictions.  One  of  the  tests  will  be  whether  the 
six  bi-monthly  correction  fields  perform  better  than  the  two 
semi-annual  correction  fields.  In  some  of  the  correction 
fields,  the  east-west  variation  is  not  large.  It  is  possi¬ 
ble  that  a  correction  dependent  only  on  latitude  may  perform 
as  well.  Such  a  correction  would  be  easier  to  apply.  The 
numerical  ocean  prediction  experiments  necessary  to  demon¬ 
strate  the  usefulness  of  the  corrections  are  in  progress, 
and  will  be  reported  separately. 
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Fig-  2  Correction  f  j.eld  (cal  cm  h  )  to  be  applied 
to'  FMOC  heat  flux  fields  based  on  all  36  of  the 
monthly  sets  of  (A  fi“Qr  1  evaluated  over  two-month  in¬ 
tervals.  Negative  values  indicate  that  the  upward 
heat  flux  is  to  be  reduced  by  the  amount  shown. 
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Appendix  A 

MONTHLY  SETS  OF  ^  H-Q^,  AND  QT 

The  monthly  sets  of  ^H-QT,  &  H  and  2r  based  on  two- 
monthly  intervals  are  included,  except  for  the  May  -  July 
1976  set  which  is  in  Fig.  1.  The  form  of  these  diagrams  is 
described  in  the  caption  of  Fig.  1.  Positive  values  (solid 
lines)  are  increasing  heat  content  in  time  and  net  surface 
heating,  and  conversely  for  negative  values  (dashed  lines). 
Zero  lines  are  enhanced. 
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